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Abstract. Bearing on the model for the time-dependent metagalactic radiation field developed in the first paper of 
this series, we compute the gamma-ray attenuation due to pair production in photon-photon scattering. Emphasis 
is on the effects of varying the star formation rate and the fraction of UV radiation assumed to escape from the 
star forming regions, the latter being important mainly for high-redshift sources. Conversely, we investigate how 
the metagalactic radiation field can be measured from the gamma-ray pair creation cutoff as a function of redshift, 
the Fazio-Stecker relation. For three observed TeV-blazars (Mkn501, Mkn421, H1426-I-428) we study the effects 
of gamma-ray attenuation on their spectra in detail. 



1. Introduction 

High-energy gamma rays traveling through intergalactic 
space can produce electron-positron pairs in collisions 
with low energy photons from the metagalactic radiation 
field (MRF) (Nikishov 1962, Goldreich & Morrison 1964, 
Gould & Schreder 1966, Jelley 1966). Gamma rays of en- 
ergy above 1 TeV typically interact with infrared photons 
of wavelengths larger than 1 /im, such as those predomi- 
nantly emitted by dust-obscured galaxies. Gamma-rays of 
energy below 1 TeV interact with near-infrared, optical, 
and ultraviolet photons, mostly from stars. 

Following the discovery of extragalactic TeV gamma 
ray sources (Punch et al. 1992), the effects of cosmological 
pair creation on their spectra have been studied by a num- 
ber of authors (e.g., Malkan & Stecker 1998, Konopelko 
et al. 1999, Primack et al. 1999). Results differ by a large 
margin, owing to the different models for the MRF em- 
ployed by the authors, and lead to extreme physical in- 
terpretations of the observed gamma-ray spectra. The de- 
bate culminates in the claim of Meyer & Protheroe (2000) 
that the weakness of the observed gamma-ray attenuation 
might have to be remedied by a threshold anomaly for the 
pair creation process, such as predicted in certain (ad hoc) 
models of quantum-gravity which violate Lorcntz invari- 
ance (Stecker & Glashow 2001). 

Discrepancies between the models for the metagalactic 
radiation field can be traced to the different formalisms 
which are used (for a detailed discussion see Hauser & 
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Dwek 2001). The simplest method (backward evolution) 
is to extrapolate present day data or template spectra 
to high redshift in a certain wavelength range (Malkan 
& Stecker 2001). Cosmic chemical evolution models self- 
consistently describe the temporal history of globally aver- 
aged properties of the Universe (Pei, Fall & Hauser 1999). 
Semi-analytical models are invoking specific hierarchical 
structure formation scenarios to predict the MRF. In our 
approach we developed (Kneiske et al. 2002; Paper I) a 
semi-empirical, forward-evolution model for the optical- 
to-ultraviolet MRF and for the infrared part a backward 
evolution model based on the data obtained from recent 
deep galaxy surveys. The model parameterizes the main 
observational uncertainties, (i) the redshift dependence of 
the cosmic star formation rate, and (ii) the fraction of UV 
radiation released from the star forming regions. 

The expected effects of gamma ray absorption vanish 
below 10 GeV (out to redshifts of z ~ 200, see Zdziarski 
& Svensson 1989). On the contrary, at energies above 300 
GeV gamma rays suffer absorption for sources at red- 
shifts z > 0.2. This is in line with current observations. 
The satellite-borne EGRET- All-Sky-Survey has resulted 
in the 3rd EGRET catalog (Hartman 1999) of 93 blazars 
with gamma ray emission between 100 MeV and 10 GeV. 
The ground-based pointing telescopes Whipple, HEGRA, 
CANGAROO, and CAT have found only 6 well estab- 
lished blazars at energies above 300 GeV searching for 
gamma ray emission from cataloged sources, in spite of 
their superior sensitivity. The few detected sources indeed 
have very low redshifts (z < 0.2). Intrinsic absorption in 
the gamma ray sources, e.g. due to the strong infrared ra- 
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diation field associated with a dust torus, could give rise 
to a redshift-independent pair creation cutoff (Donea & 
Protlieroe 2003). 

Only three blazars Mkn 421, Mkn 501, and H1426+428 
were bright enough to determine their spectra. Mkn 501 
and Mkn 421 are showing a turn-over at almost the same 
energy adopting a power-law-times-exponential model. 
The spectra show some small changes with flux, and it 
is under investigation whether this affects the cutoff en- 
ergy or not. Both sources have almost the same redshift 
z — 0.03 and should show cosmological absorption at 
the same energy of 6 — 18 TeV. In spite of large sys- 
tematic errors, several attempts to probe galaxy evolu- 
tion from the column depth of the extragalactic infrared 
photons inferred from the gamma ray data have been 
launched (Stecker & de Jager 1993, Biller et al. 1995, 
Madau & Phinney 1996, Primack et al. 1999, Renault et 
al. 2001, Kneiske et al. 2002). A third source at a red- 
shift oi z = 0.129 has very recently been discovered. The 
statistics of the signal are poor, probably resulting from 
absorption much stronger at four times the distance than 
for the other blazars (Aharonian et al. 2002, Petry et al. 
2002, Costamante et al 2003). 

If gamma-ray sources could be detected at redshifts 
z > 0.2 using imaging air Cherenkov telescopes with 
threshold energies as low as 10 GeV, such as the MAGIC 
telescope, one could infer indirectly the MRP from in- 
frared to ultraviolet wavelengths. 

The plan of this paper is to employ various parame- 
ter sets for the calculation of the MRF and to study their 
effects on gamma-ray absorption. Friedman-cosmology pa- 
rameters were fixed to the values fim = 0.3, = 0.7, and 
h = 0.65 corresponding to the ACDM cosmology. Five 
generic constellations of MRF model parameters which 
lead to MRF spectra consistent with observations, and 
which bracket the range of allowed values, are described in 
Sect. 2. In Sect. 3 we will use the optical depth of gamma- 
rays to calculate the intrinsic spectra of Mkn501, Mkn421 
and H1426-I-428 from the observed spectra. The absorp- 
tion at higher redshifts will be discussed in the fourth 
section introducing the Fazio-Stecker relation. 

2. The Metagalactic Radiation Field (MRF) 

The model developed in Paper I accounts for the emis- 
sion from stars, ISM, and dust in galaxies. Since galaxies 
strongly obscured by dust do not show up in optical galaxy 
surveys, recent infrared and sub-millimeter surveys were 
additionally taken into account. The latter surveys show 
that more than half of the cosmic star formation might be 
hidden in obscured galaxies (Chary & Elbaz 2001). 

We will briefly outline the main ideas behind the MRF 
model developed in Paper I (Kneiske et al. 2002). The 
power spectrum of the MRF in the comoving frame is 
given by 
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For the spectra L„{t) of a simple stellar population 
as a function of age r we used population synthesis mod- 
els by Bruzual & Chariot (1999) with low metalicity. The 
absorption due to interstellar medium is modelled by uni- 
form distributed dust and gas. In Paper I, the gas around 
the stars is assumed to absorb all UV photons originat- 
ing from young, massive stars. Here we introduce a new 
parameter which is the fraction of the ionizing photons 
that escapes from a galaxy, e.g. through super-bubbles 
blown into the ISM by supernovae (cf. Stecker et al. 1998). 
This fraction is, depending on the galaxy type, quite low 
(i.e. Starbursts: fesc < 6% Heckmann et al. 2001). Note 
that UV photons from active galactic nuclei (AGN), which 
might contribute considerably at UV wavelengths, are not 
considered at this stage. We adopt an average extinction 
curve for the dust absorption. The reemission is calculated 
as the sum of three modified Planck spectra 



LiiLboi) = V c,{Lboi) ■ Qx ■ Bx(Ti) 



i=l 



(3) 



where Qx oc and Li,oi = Li,oi{t), here r is the age of 
the stellar population. The three components characterize 
cold (ci, Ti) and warm (c2, T2) dust and the contribution 
due to PAH molecules (ca, T3). To get the best- fit values 
for the parameters we used a sample of galaxies detected 
with IRAS at 12/xm, 25/im, 60/im and 100/im. The warm 
dust component has its maximum around 50/Lt and is quite 
low in our best-fit model. Increasing C2, we obtain a model 
with an enhanced fraction of warm dust, the Warm-Dust 
model. 

The global star formation rate (SFR) p^:{z) consists of 
two components. 
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The first component accounts for the stars which can be 
optically detected, the other coming from stars hidden by 
dust which can only be seen looking at infrared or submm 
wavelength. Each of the SFRs can be approximated with 
a simple broken power law 



/9,(z) CX (1 + zY 



(5) 



with a = a,fn > for z < Zpeak and a = (3rn < 
for z > Zpeak- So cach SFR provides four parameters 
ttm, /3m, Zpcak and /o',(2;poak)- The values we used are 
shown in Tabled 

For reasons outlined in Paper I, it is sufficient to con- 
sider the star formation rates and the UV escape fraction 
as the dominant parameters. In the following, we will use 
six generic parameter constellations leading to MRF pre- 
dictions bracketing the data, see Fig. Q 
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Fig. 1. Comoving- frame metagalactic radiation field (including UV component) at various redshifts. "Best-fit" model, 
thick solid line; " Warm-Dust" model, thin dashed line ;"Low-IR" model, dot-dashed line; "Low-SFR" model, thin 
solid line; " Stellar- UV" model, dashed line; and " High-stellar-UV" model, dotted line. Data at z = 2, 3,4 are taken 
f 
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Fig. 2. Total global star formation rate for different mod- 
els SFR = SFRoPT + SFRl/g (for details, see Paper I) 



2.2. Warm-Dust model 

This model is almost identical to the Best-Fit model 
except for the amount of warm dust in the interstellar 
medium. The different dust types are calculated to fit the 
line intensities detected in a sample of infrared galaxies 
by fRAS at 12/im, 25/im, 60/zm and lOO^m. The range 
between 25^m and 60/im is not well determined. We raise 
the amount of warm dust at 80 K to the maximum deter- 
mined by the spectra of individual infrared galaxies. The 
corresponding change is clearly noticeable in all panels of 
Fig. n and in Fig. for sources at redshifts smaller than 
z w 0.02. 



2.3. Low-IR model 



The infrared part of the low redshift (z<0.5) MRF is 
important for the absorption of gamma rays from low 
redshift sources. In order to consider the least possible 
2.1. Best-fit model gamma-ray attenuation, we adopt a low-IR model, where 

we choose the infrared star formation rate as low as al- 
The parameters used in the best-fit model interpolate best lowed by observational lower limits on sub- mm galaxy 
the data from galaxy number counts at optical and in- number counts (e.g. SCUBA « 0.06 for 2 < z < 5 Hughes 
frared wavelengths, and direct measurements of the extra- et al. (1998) , see FigEJ- For the sake of demonstration 
galactic background, i.e. the present-day MRF. The model we accept that the present-day background in this model 
MRF shows a sharp cut-off at 0.1 /im due to the total ab- somewhat drops below the lower limits from ISO at 15/im 
sorption of ultraviolet starlight due to interstellar gas. and from IRAS at 60 /xm. 
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Table 1. Parameters (Definitions see Paper I) 



part of the MRF at wavelengths A > 0.2 fim is not exactly 
the same as in the best-fit model (although it might ap- 



p*(zp) pear so on the double- logarithmic plot washing out minute 

[Me Mpc-^ yrdeta ils) . 



Best-fit model 

SFRopT 
SFRlig 

fesc — 

Warm-Dust model 

SFRopT 
SFRlig 

fesc — 

Low-IR model 

SFRopT 
SFRlig 

fesc — 
C2 = 10-2* 

Low-SFR model 

SFRopT 
SFRlig 







C2 = 10-"" 

Stellar-UV model 

SFRopT 
SFRlig 

fesc — 1 
C2 = 10-2* 

High-Stellar-UV model 

SFRopT 
SFRlig 

fesc = 4 
C2 = 10-2* 



3.5 
4.5 



3.5 
4.5 



3.5 
4.5 



3.5 
4.5 



3.5 
4.5 



3.5 
4.5 



-1.2 



-1.2 



-1.2 



-1.2 



-1.2 



1.2 
1.0 



1.2 
1.0 



1.2 
1.0 



1.2 
1.0 



1.2 
1.0 



1.2 
1.0 



0.1 
0.1 



0.1 
0.1 



0.1 
0.03 



0.1 
0.1 



0.1 
0.1 



0.1 
0.1 



2.6. High-stellar-UV model 

The proximity effect (Bajtlik, Duncan & Ostriker, 1988) 
allows for a measurement of the metagalactic UV radiation 
field at a wavelength of 912 A, even at very high redshifts. 
Comparing our results with the recent compilation of data 
from Scott et al. (2000) (upper panels Figure ^ shows 
that even the model with fesc = 1 lies below the data. 
The reason for this discrepancy could be that the MRF 
model does not yet include the UV emission produced by 
AGN. Haardt & Madau (1996) found values for an AGN 
contribution to the UV MRF of 1.64 nW m^^ g-i sr^^ at a 
redshift of z = 2.5 which would account for the value mea- 
sured with the proximity effect method. The UV excess 
over the stellar MRF is about a factor of four and should 
make a significant difference for high redshift gamma-ray 
absorption. As a zeroth-order approximation, we assume 
an empirical UV background component by multiplying 
the UV component of the stellar-UV model by a factor 
of four, matching at the Lyman break by linear interpo- 
lation. Chosing a, perhaps more realistic, power law rep- 
resentation of the template would not make a significant 
difference for this analysis. 

2. 7. The Optical Depth 

The optical depth for pair creation for a source at redshift 
Zq, and at an observed energy Ej, is obtained from 



Tjf {Ej , Zq) 



-^^•n(2, e)-ajj{Ej, e, fi, z') de dfi dz'{6) 



2.4. Low-SFR model 

The cosmic star formation rate at high redshifts is still a 
matter of debate. Therefore, we consider a steep decline 
of the star formation rate at high redshifts (see Fig ^ as 
an extreme to study the effect of high z star formation 
and gamma ray attenuation (as opposed to the plateau in 
the other models). The star formation rate in the low-SFR 
model drops so rapidly that at z = 2 it is already an order 
of magnitude below the star formation rate of the other 
models. 



2.5. Stellar-UV model 

This model differs from the best-fit model in the escape 
fraction parameter which is set to unity allowing for the 
entire stellar UV light to escape into the metagalactic 
medium. Due to the missing reprocessed UV radiation, the 



with the cosmological line element jp- , the angle 9 between 
the interacting photons /i = cos(0), the number density of 
the MRF n{z, e) as a function of reshift and MRF photon 
energy and the pair-production cross section cr^^. 

By comparing the generic MRF models (see Fig.EJ, it 
can be seen that the optical depth from 0.2 < z < 1 is 
rather insensitive to the choice of the parameters in the 
MRF model. The interaction mainly takes place with opti- 
cal MRF photons, which are emitted by stars and undergo 
no absorption of the ISM. The rather small differences in 
pure stellar models have only a weak effect on the optical 
depth. However, at smaller redshifts effects become strong. 
For example, using the best-fit model we obtain a cut-off 
energy of ^ 5 TeV, while for the low-IR model the energy 
is ^ 18 TeV. Only at very low redshifts (z < 0.1) the effect 
of the Warm-Dust model is noticeable. The cut-off energy 
for a source dX z — 0.03 (redshifts z > 0.02) remains 
largely unaffected, see Fig. O and Fig. 01 since the change 
is relevant only for energies above 10 TeV and occurs at 
an optical depth from 2 to 10. Only Fig.Ois showing such 
high gamma energies at the high end of blazar spectra. 
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Fig. 3. Optical depth for various redshifts adopting a ACDM cosmofogy. The labeling of the line styles is the same as 
in Fig. 1. The crossing point with the line t — 1 defines the exponential cutofF energy. 



The upper limit of the shaded region above 10 TeV are 
made by the warm-dust model. Looking at spectra of sin- 
gle sources is the only method to probe this part of the 
MRF (e.g. Biller et al 1995, Stanev fePranceschini 1998, 
Guy et al. 2000, Renault et al. 2001). But there are two 
problems, the low statistic of the TeV data and the pro- 
duction processes of such high energy photons in AGNs 
which are still under discussion. 



At higher redshifts, the optical depth due to interac- 
tions with the UV part of the MRF become important. 
Consequently, the cut-off energy decreases by adding the 
UV components to the MRF model. By contrast, the cut- 
off energy increases by lowering the star formation rate at 
high redshifts in the low-SFR model. For example, the cut- 
off energy for a source at a redshift of z = 4 ranges between 
-16 GeV for the high-stellar-UV model and ~ 40 GeV 
for the low-SFR model. 



A comparison with other models is shown in Fig ^ 
The range bounded by our models contains the results of 
Dwek et al. (1998), Chary& Elbaz (2001), and Franceshini 
et al. (2001) in the infrared part of the spectrum. There are 
differences in the optical part where the other models come 
just close to the lower limits given by the HST data. At low 
redshifts, this part is irrelevant for gamma-ray absorption 
due to the threshold condition, so the difference is of no 
concern to us. 



3. Spectral modifications of nearby blazars 

A number of extragalactic gamma-ray sources have been 
detected with imaging air-Cherenkov telescopes (Table 
6, Horan et al. 2002). Three of them (with redshifts 
z — 0.03,0.03,0.129) were bright enough to resolve their 
spectra in the TeV energy band, as shown in Figure 0. 

The observed spectra are modified by gamma ray at- 
tenuation, i.e. 



(7) 



where t^^{E, z) is given by Eq. 1 (see FigureElfor a num- 
ber of examples) . Note that we neglect secondary gamma 
rays arising from cascading in the metagalactic radiation 
field, which are discussed in a separate paper (Bretz et al. 
2003). We used the best- fit model, the Warm-Dust model 
and the low-IR model to bracket the range of the unab- 
sorbed (intrinsic) spectra. Using the other generic models 
would not alter the results, since the optical depth at low 
redshifts is practically independent of the ultraviolet part 
of the MRF. 

The intrinsic spectrum of MknSOl (power per band- 
width) shows a moderate peak at 1-4 TeV. In a similar 
analysis, de Jager & Stecker (2002) found 5-9 TeV for 



1 World-data set for MknSOl: Aharonian et al. (1999, 2001), 
Kranich et al. (2001), Krennrich et al. (1999), Djannati-Atai et 
al. (1999) and Hayashida et al. (1998); for Mkn421: Kohnle et 
al. (2001), Cortina et al. (2001) and Krennrich et al. (2001); for 
H1416-I-428: Petry et al. (2002), Horan et al. (2002), Aharonian 
et al. (2002) and Costamante et al. (2003) 




Fig. 4. MRF (upper panel) and the optical depth (lower panel) at redshift z = and z = 0.03 respectively. Shaded 
regions denote the regions bounded by our models. Dwek et al. (1998), solid hne; Franceschini et al. (2001), dashed 
line; Chary&Elbaz (2001), dot-dashed line. Data references see Fig. 1. 



the peak depending on the infrared model used. They 
claim that this peak could correspond to the observed 
X-ray peak at 50-100 keV (adopting an SSC model). 
While Krennrich&Dwek (2003) found quite lower values 
for Mkn501 785 ± 153GeV to 2390 ± 127GeV. 

The intrinsic spectrum of Mkn421 still seems to show 
a shallow turnover. Assuming an exponential cutoff, the 
energy changes from 3 TeV (observed, i.e. absorbed) to 
4 TeV (intrinsic, i.e. unabsorbed). A peak is very hard 
to see, but if there is a maximum it is consistent with 
the values found by Krennrich&Dwek (2003) and would 
lie between 0.5 and 2 TeV. This could be the signature 
of the inverse Compton peak reflecting the X-ray peak 
at 6-8 keV. To further disentangle observed and intrinsic 
spectra, it is helpful to look at flux-dependent spectra, 
using the defining blazar property of being highly vari- 
able sources. The Mkn421 high-flux spectrum seems to 
be curved stronger than the low-flux spectrum which re- 
sembles a power law (Aharonian et al. 2002(1)). However, 
the statistical bias in these studies is non-negligible, and 
no conclusive evidence has emerged. A conservative esti- 
mate to study the location of the peak at different flux 
levels is introduced in Krennrich&Dwek (2003) and will 
be discussed in more detail in a paper by the same au- 
thors. They found a shift in the peak energy between the 
lowest and highest flux levels. The cut-off could also be of 



an intrinsic origin, if strong radiation fields surround the 
gamma-ray production zone (Mannheim 1993, Donea & 
Protheroe 2003). Strong radiation fields in the far-infrared 
to near-infrared wavelength ranges result from irradiated 
dust tori generally surrounding the central engine in AGN. 
However, low-redshift blazars (BL Lacertae objects) seem 
devoid of massive tori, and show sub-Eddington accre- 
tion implying weak irradiation. Moreover, it is at present 
not clear how deep in the dust torus the gamma rays are 
produced rendering this mechanism difficult to estimate. 
Theoretical SSC models (Caprini 2002) or proton blazar 
models (Muecke et al. 2003) including time variability do 
not yet predict the location of gamma ray emission zone 
relative to the dust torus. 

Inferring a spectrum for H1426-I-428 is difficult owing 
to poor statistics. Inspecting the range of the likely intrin- 
sic spectrum, the low flux seems to be the consequence of 
heavy absorption owing to the comparatively high red- 
shift of the source (four times larger than for Mkn421 and 
Mkn501, respectively). The expected exponential cut-off 
energy at a redshift of z = 0.129 obtains values of 100-200 
GeV, i.e. at energies below the threshold energy of the 
detecting instruments (Whipple, HEGRA). 

Calculating the intrinsic spectral range we find for the 
best- fit model the same result as Costamante et al. (2003). 
The energy spectrum increases with energy. However, in- 
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Fig. 5. World-data sets for three TeV blazars, and ranges of their intrinsic (" de- absorbed" ) spectra. 



spection of Figure shows that the intrinsic energy flux 
spectrum inferred from the low-IR model remains rather 
flat, or even shows a shallow downturn, implying a peak 
energy an order of magnitude lower than in the calcula- 
tion of Costamante et al. (2003). As shown in Costamante 
et al. (2001), the X-Ray peak is at an energy around or 
larger than 100 keV, and this would argue in favor of the 
gamma-ray peak larger than 12 TeV (adopting an SSC 
model). In Bretz et al. (2003), we discuss in detail the 
fate of the absorbed gamma ray photons which carry a 
substantial energy flux. 

4. The Fazio-Stecker relation 

The energy-redshift relation resulting from the cosmic 
gamma-ray photosphere T-y^{E.y,z) — 1 depends on the 
column-depth of the absorbing photons, as can be seen 
from inspection of Eq. (6). We coin this relation, plotted 
in Fig. (6), which proves to be very useful to study the 
MRF, the "Fazio-Stecker relation (FSR)" (first shown by 
Fazio & Stecker 1970)^ The theoretically predicted FSR 
(depending on the MRF model and cosmological param- 
eters) can then be compared with a measured one, by 
determining e-folding cut-off energies for a large sample 
of gamma ray sources at various redshifts. Two impor- 
tant corollaries follow from inspecting the Fazio-Stcckcr 
relation: (i) gamma-ray telescopes with thresholds much 
lower than 40 GeV are necessary to determine the cut- 
off for sources with redshifts around the maximum of star 

^ In 1968, Greisen has already suggested (in a lecture 
Brandeis Summer Institute in Physics) that pair-production 
at high redshift between optical and gamma photons would 
produce a cut-off around 10 GeV. 



formation z ~ 1.5, and (ii) gamma-ray telescopes with 
a threshold below 10 GeV have access to extragalactic 
sources of any redshift (another cosmological attenuation 
effect sets in at 2; ^ 200, Zdziarski 1989). 



The main obstacle for this method to indirectly mea- 
sure the MRF by achieving convergence between theoret- 
ical and observed FSR is the uncertainty about the true 
shape of the gamma ray spectra before cosmological ab- 
sorption has ocurred. In the simplest case, the intrinsic 
spectra would be just power law extensions of the (defini- 
tively unabsorbed) lower energy spectra to higher ener- 
gies, representative of non-thermal emission. Even in this 
optimistic case, it is difficult to assess this lower energy 
spectrum owing to the source variability. Simultaneous 
multi-wavelength observations are required. However, life 
is expected to be more complex, and the intrinsic spectra 
might consist of a sequence of humps, with spectral hard- 
ening and softening in the observed energy window (this 
is the case, for instance, in the proton-initiated cascade 
models, e.g. Mannheim 1998). Moreover, the presence of 
strong infrared radiation fields within the sources would 
lead to gamma ray cut-off energies well below the cos- 
mological ones (Donea & Protheroe 2003). These effects 
would mostly place the measured pairs of cut-off energy 
and redshift to the left of the expected FSR. The obsta- 
cles can be overcome by collecting data of a large sample 
of extragalactic gamma ray sources, and by matching the 
theoretical FSR to the upper boundary curve for all en- 
tries in the FS diagram. 
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4.1. Gamma-ray attenuation for low-redshift sources 

In the redshift interval 0.02 < z < 0.1, three of the five 
used generic MRF parameter sets show a clearly distin- 
guishable behavior (see Figure El upper panel). The three 
sets are best-fit, low-IR, and low-SFR, whereas those sets 
in which the UV fraction is varied naturally do not lead 
to measurable differences. The warm-dust model is only 
different from the best-fit model at redshifts smaller than 
0.02 (dashed-line). The published values for the cut-off en- 
ergies of the three nearby blazars seem to support the best- 
fit MRF model which entails a strong far-infrared compo- 
nent in the present-day extragalactic background. No clear 
evidence for deviations from intrinsic power law spectra 
emerges, a hardening of the intrinsic spectra would allow 
for a yet stronger FIR component in the MRF model than 
adopted in the best-fit model. However, an increase by 
more than a factor of two would be inconsistent with the 
chemical abundances in the Universe (Pei, Fall & Hauser 
1999). 

4.2. Gamma-ray attenuation for high-redshift sources 

At high redshifts, the Fazio-Stecker relation for the high- 
IR and the low-IR models converge, since the optical depth 
becomes independent of the density of infrared photons. 
Due to the large distances, which increase the column 
depth for pair production, the cut-off energies are gen- 
erally at lower energies where the threshold condition im- 
plies interactions with harder photons. Hence it follows 
that the sensitivity on the diffuse UV radiation field is 
enhanced at high redshifts. Considering UV radiation up 
to the level of the proximity data in the high-stellar-UV 
model, gamma-ray attenuation for sources at redshifts 
z > 1 increases considerably compared to the best-fit 
model, or the stellar-UV model. The low star formation 
rate at high redshifts adopted in the low-SFR model would 
reduce the UV photon density, making the universe more 
transparent for gamma-photons. However, the large drop 
of the star formation rate has been adopted only to bracket 
the range of possibilities, while latest observations rather 
indicate a plateauish behavior for the star formation rate 
at high redshifts, or even a shallow upturn. 

Note that the Fazio-Stecker relation in this red- 
shift regime is also sensitive to cosmological parameters 
(Blanch & Martinez 2001), which we have set to the cur- 
rent ACDM cosmology values. 

5. Conclusions 

Direct methods for measuring the MRF based on faint 
galaxy counts suffer from being limited to rather narrow 
wavelength ranges (introducing strong selection effects) 
and from not being sensitive to a truly diffuse compo- 
nent of the MRF. Complementary information from an 
inherently independent method, such as measuring the 
Fazio-Stecker relation is therefore of great diagnostic value 
to determine the amount of star formation occurring in 



optically-obscured infrared galaxies, to find other than 
stellar sources of UV radiation at high redshifts, and to 
determine the star formation rate at high redshifts. 

A major source of uncertainty is the shape of the unab- 
sorbed gamma ray spectra, and how much of an observed 
cut off can be attributed to cosmological absorption in the 
MRF. In the Fazio-Stecker diagram, sources with intrinsic 
gamma ray absorption would appear to lie on the left side 
of the FS relation. This bias can be overcome statistically 
by observing a large sample of gamma ray sources over a 
broad range of redshifts to find the upper bound in the FS 
diagram. To probe the present-day infrared background at 
energies larger than 10/j,m using the FS diagram, cut-off 
energies from sources at redshifts z<0.02 are needed. The 
only method to make predictions in this part of the MRF 
is to look at single-source spectra at energies larger than 
10 TeV, dealing again with all uncertainties mentioned 
above. 

Practically, this emphasizes the importance of very 
large imaging air-Cherenkov telescops such as MAGIC 
(Corinta 2001, Martinez 2003) or the proposed ECO-1000 
telescope (Martinez et al. 2003, Merck et al. 2003), which 
have the capability to discover a large number of sources, 
and which achieve low-threshold energies. 

The small range of the Fazio-Stecker relation which 
has been probed with the current generation of lACTs 
indicates consistency with the optical and infrared data 
on the MRF. Lacking detections of high-redshift sources 
with cut-off energies below 100 GeV, no conclusions about 
the UV part of the MRF can currently be drawn . 
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